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Abstract 

This paper is devoted to difFusion limits of linear Boltzmann equa- 
tions. When the equilibrium distribution function is a Maxwellian dis- 
tribution, it is well known that for an appropriate time scale, the small 
mean free path limit gives rise to a diffusion equation. In this paper, 
we consider situations in which the equilibrium distribution function is 
a heavy-tailed distribution with infinite variance. We then show that 
for an appropriate time scale, the small mean free path limit gives rise 
to a fractional diffusion equation. 
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1 Introduction 

It is well known that under appropriate scaling, the asymptotic analysis 
of collisional kinetic equations can lead to diffusion-type equations. This 
scaling corresponds to a long time scale and a small mean-free path. More 
precisely, the starting point is the following collisional kinetic equation: 

dtf + v-yj = L{f) in (0,oo) xM^ xM^, (1) 
/(0,.) = /o inM^xM^, (2) 

which models the evolution of a particle distribution function, f{t,x,v) > 
depending on the time t > 0, the position x G and a variable v G 
. This variable v usually represents the velocity of the particles, or some 
other internal degree of freedom of the particles, such as a wave vector. 
For simplicity, we take v E M^, though other spaces could be considered 
(torus for the wave vector in semiconductor or bounded set for relativistic 
particles). 

We then introduce the macroscopic variables 

x' = £x t' = 9{e)t 
and the rescaled distribution function 

nt',x\v) = f{t,x,v), 
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where e is a small parameter. The function is now solution of (we have 
skipped the primes) 

e{e)dtr + evV,r = Lir). (3) 

The object of this paper is to investigate the behavior, as e goes to zero, of 
the solutions of ([3]). It will naturally strongly depend on the properties of 
the collision operator L. 

Throughout this paper, we will assume that L is a linear Boltzmann 
operator (sometimes also called "scattering" operator), describing the inter- 
actions of the particles with the surrounding medium, of the form: 

L{f) = [ [a{v, v')fiv') - aiv', v)f{v)] dv' (4) 

with a non-negative collision kernel a = a{v,v') > 0. The operator L is 
conservative, i.e. it preserves the total mass of the distribution. Under 
classical assumptions on the collision kernel a, there exists a unique positive 
normalized equilibrium function F: 

F = F{v) > a.e. on M^, / F{v) dv = I and L(F) = 0. 

In the sequel, we always assume that F exists and is an even function of v. 

The derivation of diffusion-type equations from kinetic equations such as 
([3]) was first investigated by E. Wigner [20], A. Bensoussan, J.L. Lions and G. 
Papanicolaou in ^ and E.W. Larsen and J.B. Keller |14]. In [7|, P. Degond, 
T. Goudon and F. Poupaud consider very general collision operators of the 
form ([4]). When F decreases "quickly enough" for large values of \v\ (and 
under a few additional assumptions on a and F), they prove in particular 
that for 6{e) = e^, f^{t,x,v) converges, when e goes to zero, to a function 
of the form p{t, x)F{v) where the density p{t, x) solves a diffusion equation 

dtp - {D V^p) = 0, (5) 

with diffusion matrix D given by the following formula 

D = I (v^x) dv with L{x) = -vF. (6) 

In order to introduce our problem in simple terms, we now consider, 
in this introduction and in the next section, that L has the following very 
simple form: 

LU)=pF-f, p={f):= [ f{v)dv, (7) 
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which corresponds to the choice of a colhsion kernel a{v,v') = F{v). This 
baby model is usually called "linear relaxation" or sometimes "linear BGK" 
collision operator. We shall come back to more general (and realistic) col- 
lision kernel a{v,v') in Section [3l Under d?]), it is readily seen that we can 
take X = V F in the definition of the diffusion matrix D and thus: 




In particular, in order for D to be finite, we need the second moment of F 
to be finite: 

/ \v\'^F{v)dv < oo. (8) 
This is the case, for example, when the equilibrium F is given by the so- 

2 

called "Maxwellian" distribution F{v) = Cexp(— ^). The main result of 
this paper can then be summarized as follows: If F is such that ^ does not 
hold, then, under appropriate time scale, the limit e — > in ([3D leads to a 
fractional diffusion equation instead of ([5]). 

More precisely, we will assume that F is a heavy-tailed distribution func- 
tions, that is typically (we will consider slightly more general F later on): 

^(^) ~ \v\N+a 1^1 ^ ^' 

with kq > and a > (this last condition guarantees that F is integrable). 

Heavy-tailed distribution functions arise in many contexts. For instance, 
most astrophysical plasmas are observed to have velocity distribution func- 
tions exhibiting power law tails (see Summers and Thorne [18] or Mendis 
and Rosenberg [IGj). Dissipative collision mechanisms in granular gases can 
also produce power law tails: see for instance [TTj for the so-called "inelastic 
Maxwell model" introduced in [5]- One can also refers to the more general 
review paper |19j . We also mention that a recent work [5] has shown that 
even elastic collision mechanisms can produce power law tail behaviors in 
the case of mixture of gases with Maxwellian collision kernel. Power law 
tails are also common in economy where they are referred to as Pareto dis- 
tributions (as well as in statistics and probability more generally). Here 
are some samples of (very different) mathematical works using statistical 
physics models to account for these power laws in economy: Newman [17J, 
Duering-Toscani [lOj and Wright [21]. Let us just mention to the reader 
that, in the case of economy, the interpretation of these power law distri- 
butions as equilibria of statistical physics models (which is the core of the 
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field now known as "econophysics" ) is controversial since it does not take 
into account any individual rationality. 

It is thus the goal of this paper to investigate what happens to diffusion 
limits when the velocity repartition is no longer described by a Gaussian 
function. More precisely we shall be particularly interested in those heavy- 
tailed distributions whose variance is infinite, since in this case we shall 
show that the diffusion limit yields equations of the form (jlip . and our 
underlining motivation was indeed to provide a microscopic derivation for 
these fractional diffusion equations, which is lacking at now. 

As a striking link between these goal and motivation, let us mention that 
a famous case of such distribution with infinite variance are the stable (or 
Levy) distributions. And the latter plays an important role in probability 
theory since it can be interpreted as the law of a "Levy walk" whose law 
evolution is governed by a fractional diffusion equation. To say it differently 
and in a more analytical way, these stable distributions are the fundamental 
solutions of the fractional diffusion equation. They thus play the role played 
by the gaussian distribution in the case of the heat equation. 

When ([9]) holds with a > 2 then ([8]) is still satisfied, and the analysis 
leading to ([5]) can be performed. We are thus interested in values of a less 
than 2. More precisely, when a £ (0,2), we will prove that the appropriate 
diffusion scaling is given by 

0(e) :=£°, (10) 

and that the solution of ([3]) then converges to p{t, x)F{v) with p solution 
of the following fractional diffusion equation 

+ K (-A^)"/^ = in(0,oo)xM^, (11) 
p(0,.)=po inM^, (12) 

with ^ 

— — ^ 1 — r^, — dw. (13) 

We recall that the operator (— A^;)"/^ denotes the fractional operator, de- 
fined for instance by the Fourier formula 

{-/\^r/^p:= T-\\k\^ Hp){k)), 

where J- stands for the Fourier transform in the space variable. This opera- 
tor has a lot of nice properties, similar to that of the usual Laplace operator, 
most notably it retains some ellipticity (in the sense that its Dirichlet 
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form controls some fractional Sobolev norm H"^'^). The fundamental differ- 
ence is the fact that this operator is nonlocal for any < a < 2. 

Before we state our first result, let us introduce some functional spaces 
definitions. We denote by I/, p € [l,+oo] the usual Lebesgue space on 

or X (this will always be clear from the context), with the flat 
measure. When B > is a given positive locally locally integrable weight 
function, we write U'(Q) for the Lebesgue space with measure Qdx, that is 
the Banach space defined by the norm 

ll/llL.(e):=f / I/re 

We shall sometimes use subscript letters in order to recall and emphasize 
which variable we are considering for the Lebesgue space. 

Our first result could thus read as follows: 

Theorem 1.1. Assume that L is given by ^ and that F satisfies (0) with 
a G (0, 2). Assume furthermore that fo € L?'{F~^) and let be the solution 
of ^ with 9{e) = and initial condition /q. 

Then, when e goes to zero, f^ converges in L°°(0, T; L^(M^ xMJ^))-weak 
to pF with p = p{t, x) the unique solution to the fractional diffusion equation 

m, m, m- 

Note that under ([9]), a can be characterized as follows: 

a = supia>0; / {vl^Fiv) dv < oo\ . (14) 
I Jm J 

This characterization will be more obvious in Section [3] when we add to F 
a slowly varying function. In particular, we will see that the fact that the 
moment of order a is finite or not does not affect the asymptotic equation. 

It is also worth pointing out that when L is given by ([7]), fractional dif- 
fusion is only observed when the moment of order 2 is unbounded, i.e. when 
the energy associated to the equilibrium distribution function F is infinite. 
We will see, however, in Theorem 13.21 that fractional diffusion may arise 
even when the energy of F (and higher order moments) is finite, provided 
the collision frequency 

^{v) = / cr{v' , v) dv' 

in dH) is degenerate for high velocities: u{v) ~ \v\^ for large v with /3 < 
small enough (for instance with a > 2 and /? < 2 — a < 0). 
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Conversely, even when the energy of the equihbrium if infinite, it is 
possible in some cases to recover a classical diffusion limit. This is the 
case for instance when a € (1,2) and 2 — a < (5 < 1, or a G (0,1) and 
1 < (5 < 2 — a (see Theorem 13. 6p . These ranges of the parameters /3 
are quite surprising as well as the fact that 7(0, (5) (the order of the limit 
fractional diffusion operator) is decreasing (to 1) as a function of (3 when 
a G (0, 1) and increasing (to 1) as a function of (3 when a G (1,2), whereas 
one would expect at first guess from physics that the stronger /?, the stronger 
the diffusion at the limit as it enhances collisions for high velocities. This 
calls for a satisfying physical interpretation. 

Theorem 11.11 states that in the long time limit, the particles evolve ac- 
cording to an anomalous diffusion process (one calls a diffusion process 
anomalous if the mean square displacement grows like a nonlinear func- 
tion of the time variable in the long time limit). Anomalous diffusion limits 
for kinetic models are well known in the case of a gas confined between two 
plates, when the distance between the plates goes to (see [6], [l2], [8], [9]). 
In that case, the limiting equation is still a standard diffusion equation, but 
the time scale is anomalous (0(e) ~ e^ln(e~^)). The particles travelling in 
directions nearly parallel to the plates are responsible for the anomalous 
scaling. We will see that a similar behavior arises here when a = 2 in ([9]) 
(i.e. when the second moment is unbounded, but all moments of smaller 
orders are bounded), see Theorem 13.61 

A fractional diffusion equation has been obtained as a diffusive limit 
from a linear phonon-Boltzmann equation simultaneously and independently 
by M. Jara, T. Komorowski and S. 011a in [13], by a different probability 
approach. Let us also mention a work in progress [H] by B. Texier and the 
first author, where a derivation of fractional diffusion equations from kinetic 
models is obtained in the framework of a gas confined between two plates 
with singular equilibrium distribution functions, by analytic means. 

Our approach relies on the use of Laplace- Fourier transforms and a care- 
ful computation of the asymptotic behavior of the symbol of the differential 
operator. This approach is quite simple and allows for an explicit compu- 
tation of the coefficients of the asymptotic equation under optimal assump- 
tions. Preliminary computations seem to indicate that it can generalize to 
equations involving an external force field, but it seems harder to generalize 
to a non-linear collision kernel (Pauli's statistic for instance). A more "non- 
linear" approach, as was the moments method in the classical diffusion case, 
would be welcomed. 

Let us now briefly outline the contents of the paper: The proof of The- 
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orem ll.H which corresponds to the simplest case, is presented in the next 
section. In Section 3, we state our main result (Theorem l3.2p which general- 
izes Theorem II. II to a large class of linear collision operator L of the form 
We also address the critical case (a = 2) (Theorem l3.4p and the classical dif- 
fusion case in Theorem 13.61 Theorem 13.21 (as well as Theorem 13. 6|) are then 
proven in Section H] while Section [5] is devoted to the proof of Theorem 13.41 

2 Proof of Theorem 11.11 (simplest scenario) 

The goal of the section is to present the main ideas of the proof, which are 
in fact quite simple. We thus make all sort of simplifying assumptions in 
order to focus on the important aspect of the proofs. We recall that L is 
given by ([7]), and we assume, instead of ([S]), that the equilibrium function 
F satisfies: 



F{v) < Ko |?;|-^-" for all v G M^, F{v) = kq |?;|-^-" if \v\ > 1, (15) 
with kq > 0. 



The proof is divided in 4 steps. Note that the proof of the main result 
in Section [3] will follow exactly the same steps. 

Step 1: A priori estimates for 



The solution f of ([3]) (with 9{e) 



e") satisfies 




from which we deduce the two estimates 




(16) 



and 



Jo JR'^^ 

Cauchy-Schwarz inequality also gives: 




(17) 
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so that p^{t,x), as well as L{f^), are well defined a.e., and 

sup / p%t, .f dx < Wfoh^^F-^y (18) 

Step 2: Another formulation of the rescaled equation. 

We denote by the Laplace-Fourier transform of with respect to t and 

X, defined by 

'Pip, k,v)= / e-i"' e"'^-^ f{t, X, v) dt dx, p > 0, k e M^. 
Jq 

The function then satisfies 

e"pP-e''h + eivkp = (P)F-p, 

where /o denotes the Fourier transform of the initial datum /q. We can 
rewrite this equality as 

/ = IT-;; — r— — r p + 



l + e°p + £i?;-fc l + e"p + eii'-/c' 

with p^{p, k) = {f^){p, k) the Laplace-Fourier transform of p^ , and integrat- 
ing this equality with respect to v, we obtain: 

VAjv 1 + £"^? + £^u • A; / \^7]RAr 1 £° p e z u • ^ 

The normalization condition for F now yields: 

/o 



/ 



with 



dv + (f = 0, (19) 



a'{p, k):= — [ (- : - 1^ F{v) dv. 

Let us prove that that the first term converges to po when e goes to zero. The 
assumption /o G L^(F~^) implies in particular that /o G L^(Lj). Hence its 
the Fourier transform /o also belongs to L|(LJ) by Parseval equality, which 
means that /o is intcgrablc in v for almost all k. This allows to apply the 
Lebesgue dominated convergence theorem, which yields, for almost every k, 



[ ^° ■ , dv^ [ /o 

Jt^n 1 + £'^p + £IV ■ k JkAT 



dv = Po- 



9 



So we are left with the task of studying the hmit, as e goes to zero, of 
the coefficient a^. 

Step 3: The cornerstone argument of the proof: where fractional 
diffusion symbol appears. 

A simple computation leads to 

r\ 1 /' e'^p + eiv-k 

a%p,k) = — - / -F{v)dv 

e" J^N 1 + e'^p + eiv ■ k 

f 1 + £> 

= ~P / -, —rrrivjav 

Hrn {I + e^pf + e^{v ■ kf ^ ' 

' ^ ' F{v)dv 



e° J^N (1 + e"p)2 + e'^{v- kf 

(the term involving eiv ■ k on the numerator vanishes thanks to the symmetry 
of F). The first term in the right hand side is bounded by \p\ (uniformly in 
e) and the dominated convergence theorem readily implies that it converges 
to —p Jjgjv F{v) dv = —p as e goes to zero. So it only remains to study 

If the second moment of F is bounded, then it is readily seen that with 
a = 2, we have 

d%p,k) — > {vkf F{v)dv = K\kf, kG(0,oo), 

and the limit e ^ in ([T|) leads to the diffusion equation ([5|). 

With ([15]) and a E (0,2), the second moment of F is unbounded. We 
then claim that for any p > 0, A; G M^, we have: 



\d^{p,k)\<K\k\'' and d'{p,k) — >K\k\'', (20) 

e— >0 



with K G (0, oo) given by (fT3]) . 

As a matter of fact, the first inequality in ([20|) follows from (fT5|) : 

JlRiv 1 + e^(T; • fc)^ 

where the last equality is obtained by making the change of variables w :- 
e\k\v. And in order to get the convergence of d^ , we simply write 

d'{p,k) = d\{p,k) + dl{p,k) 
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with 



< e^-" / (i; • ky F{v) dv 

J\v\<l 

< e2-«|A:p [ F{v)dv 

J\v\<l 



and 



= L>e \H (1 + + wl ^ ^ '^1"' 

where we make again the change of variables w := £\k\v and use the domi- 
nated convergence theorem. We have thus shown: 

Proposition 2.1. //a € (0,2), then 

a^{p, k) — > —p — Klk]"^ as e ^ 
with K G (0, oo) given by Furthermore, a^{p,k) satisfies 

\a'{p,k)\ < IpI + 

Step 4: Conclusion. 

From the bound (jlSp and up to extraction of a subsequence, we know that 
there exists G L°°(0, oo; L^(M^)) such that ^ rj weakly in L°°(0,oo; 
L^(M^)). On the one hand, from (jlSp again, we have bounded in 
L°°(a, oo; L^(IR^)) for any a > 0. On the other hand, for any (p = ip{t) G 
P(0, oo) its Laplace transform dp belongs to L^(0, oo) and for any ip = 
■ip{x) G SS{M.^) its Fourier transform Tip belongs to SS{M.'^) so that 

{p^, p®tp) = {p'^^C-p ® T^)) {t], Cp (g) T^) = {r],p>^ ip) 

as e — > 0. We easily deduce that p^ ^ fj weakly in L°°(a, oo; L^(R^)) for 
any a > 0. Gathering all the convergence results established above, we may 
pass to the limit in (fT9]) and we get that rj satisfies 



Po + {—p — K'\k\°') rj = for a.e. p > 0, A; G 
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Finally the unique solution p to the fractional diffusion equation pT]) , ([12]) , 
([13]) also satisfies 

Po + «:|^r)p = for a.e. p > 0, G M^, (21) 

so that rj = p a.e., and then rj = p because the Laplace- Fourier transform is 
a one-to-one mapping (say in S'{[0, oo) x M^)). We conclude that ^ pF 
weakly in L°°(0, T; L2(mA^ x M^)) thanks to ([TT]). □ 



3 Diffusion limit results for general collision oper- 
ator 

Having exposed the main ideas of the proof in the previous section, we now 
come back to more general collision operators, and state our main result. 
We recall that L now is of the form 

' dv' 



L{f) = / a{v,v')f{v')-a{v',v)f{v) 
= K{f)-vf 

with 

K{f)= [ a{v,v')f{v')dv', i^iv)= [ a{v',v)dv'. 

Such an operator is obviously linear, well-defined as a (possibly) unbounded 
operator with domain L^{i'), and closed (in fact it is bounded from L^{i^) 
to L^). It is also straightforwardly conservative: 



/ L{h) dv = for all h G L^{u). 

Jrn 



The choice of the cross-section a{v, v') is crucial. We start with the following 
structural assumption. 

Assumptions (Al) The cross-section a is locally integrahle on M^-^, non 
negative and the collision frequency v is locally integrahle on and satisfies 

v{-v) = v{v) > for all v G . 

Assumptions (A2) There exists a function < F £ L^{i^) such that 
|f P iy{v)~^ F is locally integrahle and 

v{v)F{v)=K{F){v)= I a{v,v')F{v')dv', (22) 
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which means that F is a equilibrium distribution (i.e. L{F) = 0). Further- 
more, the function F is symmetric, positive and normalized to 1: 

F{-v) = F{v) > for all v G and / F{v) dv = 1. 

Jrn 

Note that under classical assumptions on a, the existence of an equilib- 
rium function is in fact a consequence of Krein-Rutman's theorem (see [7] 
for details). A particular case in which this condition is satisfied is when a 
is such that 

yv,v' e a{v, v') = b{v, v')F{v), b{v' , v) = b{v, v'), (23) 

for some b G ^/^^(M^^). In that case, we say that a satisfies a detailed 
balanced principle or a micro-reversibility principle, while the more general 
assumption (j22p is called a general balanced principle. 

Next, we need to make precise the behavior of F and u for large \v\. 
For that purpose, we recall that a slowly varying function is a measurable 
function i : M_|_ M such that 

£{Xs) ~ £{s) as s ^ oo for ah A > 0. 

Example of slowly varying functions are positive constants, functions that 
converge to positive constants, logarithms and iterated logarithms. 

In our main result, we assume that F is a regularly varying function 
of index —{N + a) with a > and that u behaves in the large velocity 
asymptotic like a power function. More precisely, we make the following 
assumption: 

Assumptions (Bl) There exists a > and a slowly varying function I 
such that 

F{v) = Fo[v)l{\v\), (24) 

where Fq is such that 

\v\°'^^Fo{v) — > Ko G (0,oo) as \v\ oo. (25) 

Assumptions (B2) There exists /? G M and a positive constant uq such 
that 

\v\~^v{v) — s- vq as \v\ oo. (26) 
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Assumptions (B3) Finally we assume that there exists a constant M such 
that 

/ F'^dv' +( —^dv'] <M for all v € M^, (27) 
Jrn b \J^N v' ) 

with b = b{v, v') := a{v, v') F~^{v). 
Remarks & Examples 3.1. 

(i) The constant kq in i25\) could actually be a function of the direction v/\v\ 
without any additional difficulties. We will take kq constant in order to keep 
things simple. 

(a) The condition (B2) implies that the collision frequency satisfies: 

vi {vf < i^{v) < V2 {vf for large v € M^. (28) 
for some constants i^i, 1^2 G (0, 00). 

(Hi) The conditions (B2) and (B3) are fulfilled for a collision kernel a satis- 
fying the detailed balance principle 123\) where b satisfies 



b{v,v') = {vf {v'f or b{v,v') = {v-v'f, f3 < 



a. 



since then ^{v) ~ {v)^ in both cases (for the second example, we refer to the 
proof of Lemma \6.1\ in the appendix where the main arguments of the proof 
of that statement is presented). Let us notice that such collision kernels 
satisfy ^ Assumptions (Al), (A2), (A3)] but of course not in general Q 
Assumptions (B2), (B3)] except when (3 > 2 — a. 

(iv) The conditions (B2) and (B3) are also fulfilled for a (more physical) 
collision kernel a satisfying the detailed balance principle ^23\) with 

b{v,v') = \v-v'f 

and under the additional restriction (3 € (— min{a; A^/2}, min{a; A^)}. We 
refer to Lemma \ 6.1\ in the appendix where the proof of that statement is 
presented. 

(v) Our assumption (B3) is a bit more general than the corresponding as- 
sumptions (A3) in since for instance it is fulfilled by the collision kernels 
of the point (iv) above with P G (0, min{a; A^}), while such a collision kernel 
does not satisfy Q Assumption (A3)]. 

We can now state our main theorem: 
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Theorem 3.2 (Fractional diffusion limit). Assume that Assumptions (Al- 
A2) and (B1-B2-B3) hold with a > and (5 < min{a; 2 - a}. Define 



7:=^—^, and e{e) := l{e'^) e\ 

It is worth noticing that we have here P < 1 as well 7 < 2 for these ranges 
of the parameters a and (3. Assume furthermore that fo € L'^{F^^) and let 
f^ be the solution of (0^, with that choice of 9 and initial data /q. 

Then, (/^) converges in L~(0, T; L'^{R^ xR^))-weak to a function p{t, x) F{v) 
where p{t, x) is the unique solution of the fractional diffusion equation of or- 
der 7; 

dtp + K{-A^y/'^p = m(0,oo)xM^, (29) 
p(0, .) = po in R^, (30) 



with K given by 



Remarks 3.3. 



dw. (31) 



(i) Note that in Theorem \3.2\ we always have f3 < 1 from the assumptions, 
and the condition (3 < minja; 2 — a} is equivalent to the condition 7 < 2. 

(a) When (3 = (i.e. v bounded below by a positive constant), then we have 
7 = a, and we recover a result similar to that of Theorem with more 
general collision operator and equilibrium states. Note in particular that in 
that case, the addition of a slowly varying part in F has modified the time 
scale 9{e), but not the limiting equation. 

(Hi) When a G (0,1), the function (3 £ {—00, a) ^ l{l3) = {a — (3)/{l — 
(3) is decreasing from 1 to 0. When a = 1, it is defined on (—00, 1) and 
identically constant to 1. When a > 1, the function j3 G (—00, 2 — a) ^ 
7(/5) = (o^ ~ /5)/(l — /3) is increasing from 1 to 2. It is thus always possible 
to obtain fractional diffusion limit for a kinetic collisional equation with 
regularly varying equilibrium F of index — {N + a). In particular, fractional 
diffusion behavior can arise even when F has finite energy (a > 2). 

The case 7 = 2 (which may occur when a > 1) is critical in the sense 
that even though the second moment may be infinite (and the usual method 
yield an infinite diffusion coefficient), the asymptotic behavior can still be 
described by a standard diffusion equation under the appropriate time scale. 
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In this case, the exact behavior of the slowly varying function i is crucial, 
as it may determine whether the second moment is finite or infinite. When 
it is finite, (for instance £{\v\) = {ln\v\)~'^), then the usual technic yields 
the diffusion equation ([5]) under the classical time scale 6{e) = . When 
it is infinite, then the asymptotic behavior is still described by a diffusion 
equation, but the time scale has to be modified. These two situations are 
included in the next theorems: the first one corresponds to anomalous diffu- 
sive time scales, while the second one corresponds to classical diffusive time 
scales. 

Theorem 3.4 (Classical diffusion limit with anomalous time scale). Assume 
that Assumptions (A1-A2) and (B1-B2-B3) hold with 

a > 1 and f3 = 2 — a {i.e. 7 = 2), 

and i such that 

£{r) ln(r) — > -|-cxd as r ^ +00 (32) 

(note that this implies in particular that the second moment of F is infinite). 
Then define 

e{e) = e^i{e'^) \n{e-^). 

Assume furthermore that /o G L?'{F~^) and let f^ be the solution of ^3^, 
with 9{e) defined as above and initial data fo- 

Then, {f^) converges in L°°(0, T; L^(R^ x M.^))-weak to pF where p = 
p{t, x) is the unique solution to the standard diffusion equation 

dtp - K A^p = 

with K given by 

_ KqUq 1 f wj 1 

" (1-/3) n In(A-i) y|^|>, ui + wf \w\^+^ ^"'^ 
Remarks 3.5. 

(i) This is a critical case, in which the variance of the equilibrium distribu- 
tion F ( and therefore the classical diffusion coefficient ^ ) are infinite, but 
the asymptotic symbol still is of order 2. This case is still referred to as 
anomalous diffusion, even though we recover a standard diffusion equation, 
because the time scale is not the usual time scale. 
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(ii) Proceeding as Sectionl^ with a = 2, (3 = 0, and i = 1, we immediately 
see that if we would take a time scale 9{e) = , then we should expect to 
find for the diffusion coefficient: 

aw = oo. 



w 



N+2 



This is why a different time scale has to be considered in this limiting case 
a = 2. 

Theorem 3.6 (Classical diffusion limit with classical time scale). Assume 
that Assumptions (A1-A2) hold as well as the following bounds 



i'{v) \v' 



/l2 



+ ^-^ ) F' dv' <M Vug M^. (33) 



b{v,v') u{v') 

Assume furthermore that /q G LP'{F^^) and let f^ be the solution of ^3^, 

with 6{e) = and initial data /q. 

Then (/^) converges in L°°(0, T; L2(]R^ x M.'^))-weak and in L'^{{0,T) x 
X M.^))-strong to a function pF where p = p{t,x) is the unique solution 

of the standard diffusion equation ^ explained in the introduction ( with the 

same constant). 

Remarks 3.7. 

(i) Remark that Asumptions (B1-B2-B3) with any (3 > 2 — a (i.e. 7 > 2j 
obviously implies liS^) . so that all the values of a > 0, P < min{l;a} are 
addressed by Theorems \3.2l \3.4\ and [3^ One also sees from these exemples 
of coefficients that this theorem covers cases where the velocity distribution 
has infinite variance. 

(ii) Theorem \3.6\ is (in some direction!) slightly more general than ^ The- 
orem 1] since for instance /?, Asumption (A3)] implies that b is uniformly 
bounded from below, whereas Theorem \3.6\ applies to the collision kernels of 
Remarks & Examples \3.1\ (iv) with (3 G (max{2 — a; 0}, min{a; A^}), while 
they are not dealt with by [1, Theorem 1]. 

Figure 1 summarizes a large part of the results presented in the last three 
theorems. 

We now turn to the proofs of Theorems 13.21 and 13.41 presented respec- 
tively in Sections S] and [5j Concerning Theorem 13.61 we shall not write a full 
proof since it just follows from the proof of [TJ Theorem 1] by replacing the 
proof of the auxiliary result in ^ Proposition 1] by Lemma |4. II below. 
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Figure 1: Summary of the main results 
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4 Proof of Theorem D 



The general idea of the proof is the same as that of Theorem 11.11 (see Sec- 
tion [2]) . The main difference is in the derivation of formulation (jl9p : It is 
no longer possible to work with the density directly, and we use K{f^) 
instead leading to the corresponding formulation (j39p . This leads to another 
difficulty, since (j39p involves an additional term which we have to show is of 
smaller order. Finally, we will see that the computation of the asymptotic 
symbol is a little bit more complicated than in Section [2] because of the 
collision frequency being velocity dependent, and because of the presence of 
the slowly varying function t m. F (this is the object of Proposition 14.51 in 
which we see appearing the importance of the different time scale that takes 
into account i and (3). 

Throughout this section (and the next one), we denote 

(/.(e) :=£(e-T^). 



4.1 A priori estimates 



The next lemma summarizes the key properties of the collision operator L: 
Lemma 4.1. The operator is bounded in L'^{uF~^) and satisfies: 

I L{f)ldv<-^[ \f-{f)F\'^dv for all f L\vF~^) {M) 

where (/) = J^m f{v) dv and M is defined in (27^ . 

Proof of Lemma 14.11 We adapt the proof of [TJ Proposition 1 & 2]. 

To show that is bounded in L'^{h'F~^), we obviously only have to 
check that -K is bounded. Using Cauchy-Schwartz inequality and the fact 
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that K(F) = uF, we get: 



-Kif) 

V 



< 



, M v)F{v) I a{v,v')^-^f^dv' dv 



In order to prove (|34|) . we write 

L{f)^dv = [ [ a{v,v')f'^dvdv'-[ v{v)^ dv 

f f f P 

a{v,v')F' —— dv dv' — / v{v)—dv. 

F F Jm F 

Next, we note that the second term in the right hand side can be rewritten 

p2 r r f2 




J 

u{v)ydv 



[ [ a{v',v)F^dvdv' 

[ [ cr{v,v')F'^dvdv, 
Jrn J^n F'^ 



as well as (using the fact that uF = K(F)) 

f2 r f2 



f u{v)Ldv = [ K{F)Ldv 
Jri^ F Jrn F^ 

= 11 aiv,v')F'^dvdv. 



We deduce 

L{f)Uv = -\f [ <y{v,v')F' 
^ ^ Jrn Jrn 

Integrating (in the v' variable) the identity 



F' F 



dvdv'. (35) 



fF'-fF=(^-^) FF' 
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we get 



where g = f — {f)F- The Cauchy-Schwarz inequahty imphes 



so that 



/ f F ..\ f f ff /'A' WW 



^udv<{supu -F'dv'\{ / [^-^] aF'dvdv 



(36) 



Gathering ([27]), ([35]) and ([36]) we obtain □ 

Using Lemma l4.ll we can prove the following estimate on 

Lemma 4.2. The solution of ^ is bounded in L°°(0, oo; ^^(F^^)) uni- 
formly with respect to e. Furthermore, it satisfies: 

f = p'F{v)+g\ 

where the density = J^pf f^ dv and the function are such that 

IIP^IIl°°{0,oo,L2) < II/o||l2(f-i) (37) 

and 

ll5'^llL2(0,oo;L2(i/F-i)) < II /o ||l2(F-1) ^(s)^''^ • (38) 

In particular converges L°°(0, T; L^)-weak to p, and converges L°°(0, T; L?'{F 
weak to f = p{t,x)F{v). 

Proof of Lemma 14. 2L Multiplying ([3]) by /F, we get: 

/■ \f^\^\,dxdv = ^ [ L{f')^ 
2 dt L2N '-^ ' F e e) L2N ' F 



2Me{. 

which gives: 



\f^(t,x,v)\—dxdv + 



If 1 

2 ./in>2JV r 
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This inequality shows that is bounded in L°°(0, oo, L^(F ^)). Further- 
more, denoting = — p^F we also get 

/ / l/p^dxdvds < C||/o||l2(^-i) 61(e). 

Jo iM2JV i' 



Finally, Cauchy-Schwarz inequality implies: 




4.2 Another formulation of the rescaled equation 

Proceeding as in Section [21 we denote by f''{p,k,v) the Laplace- Fourier 
transform of f^{t,x,v) with respect to t and x, defined by 

f^{p,k,v)= / e-P^e-'^-''fit,x,v)dtdx. 
Jr'^ Jo 

We define p^{p,k) and g^{p,k,v) similarly, and we denote by fQ{k,v) the 
Fourier transform of /o with respect to x. 

Now, taking the Laplace- Fourier transform in ([3]), it is readily seen that 
satisfies 

e{e)pj^ -e{s)fo + iv-kT' = K(p)-'^T', 
which easily yields 

) + 0{e)p + eiv ■ k u[v) + d{£)p + eiv ■ k 

Multiplying this equality by a{w,v) and integrating with respect to v, 
we obtain 

i^f7^\f \ f 0{e)a{w,v) 7r 
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Finally, integrating with respect to w, we get: 

e{eMv) 



K{f<'){v)dv 



and thus 
: 



+ 



u{v) 



^{v) + 9{e)p + eiv ■ k 



v{v) + 9{e)p + eiv ■ k 



i^{v) + 9{e)p + eiv ■ k 
fo{v) dv 



fo{v) dv 

K(p){v)dv, 



+ 



1 



eie) 



v{v) 



v{v) + 9{e)p + eiv ■ k 



K{f^){v)dv 



Next, using Lemma |4.2| we write 

which leads to (using the fact that K{F) = vF): 

K(p)=?yF + K{g^). 

We deduce: 

v{v) 







/ 



h'{v) + 9{e)p + eiv ■ k 



fodv 



1 



9{e) 
1 



9{e) 



v{v) + 9{e)p + eiv ■ k 



v{v) + 9{e)p + eiv ■ k 



uF dv I 



1 



K{(f) dv 



(39) 



The rest of the proof consists in passing to the limit e ^ in (|39p . 
In the next two subsections, we shall show that the last term vanishes 
in Lemma 14.51 and we study in Lemma 14.31 the limit of the second term 



a^{p,k) :-- 



1 



9{e) 



i^{v) 



v{v) + 9{e)p + eiv ■ k 



1 



V F dv 



as e goes to zero. This limit will provide the Fourier-Laplace symbol of the 
asymptotic equation, and is therefore the cornerstone of the proof. 
Let us prove that the first term converges to po 



u{v) 



V + 9{e)p + eiv ■ k 



fodv 



Po, 



for almost all k G 



p G 



(40) 
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The assumption Jq G L'^{F ^) implies in particular that /o G L'].{L\). 
Hence its the Fourier transform /o also belongs to L|(L^) by Parseval equal- 
ity, which means that /o is integrable in v for almost all k. Together with 
(j4ip . this allows to apply the Lebesgue dominated convergence theorem, 
which yields, for almost every k, 



v{v) 



V + 9{e)p + eiv ■ k 



fodv 



fo dv = po 



which proves (jlO 



4.3 Passing to the limit in a"^ 

The main goal of this section is the proof of the following proposition: 
Lemma 4.3. Recall that 7 < 2 and 9{e) = e^ip{e), and let 



Jrn 



V + 0(e) p + eiv -k 



1 



u(v)F(v) dv. 



Then 



a^{p, k) — > —p — K\k\^ 

when £ goes to zero, locally uniformly with respect to p > and k G M.^ , 
with K given by I131\) . Furthermore, is locally bounded in [0, 00] x M.^ 
uniformly w.r.t. e: There exists a constant C such that 

\a%pM < \p\ + C(l + \k\^). 



Proof of Lemma 14.31 Observing that 



< 2, 



V + 6{£)p + eiv ■ k 
we see that is well defined for any e > 0. Next, we write 

v + e{e)p 



1 



u + 6{e)p + eiv ■ k {v + 9{e)p)'^ + [sv ■ k)"^ 

e^vkf 



e{e)p 



+ 
+ 



{v + 9{e)pY + {ev ■ ky 

euiv ■ k 
{y + e{e)pf + {evkf' 



(41) 



(42) 
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Using the fact that F{—v) = F{v) and I'i—v) = v{v), we deduce 



^ ^ (e^^-A:)^ dz;. (43) 



9{e)J^M {u{v)+e{e)pY + {evkf 

The dominated convergence theorem immediately impUes that the first term 
in the right hand side converges to —p F{v) dv = —p as e goes to zero. 
Furthermore, that term is clearly bounded (in absolute value) by \p\. So it 
only remains to show that (and here we recall that 6(e) := £^ip{e)) 

converges to k|A;|''' and is locally bounded when e goes to zero. 
For some M > 0, we write 

d'ip,k)=dl{p,k)+dUp,k) 

where 

dl{p,k) = [ —— ^^"""^^^ , u{v)F{v)dv 

1^^' ' e-i^ie) J\,\<MHv)+e^'f{e)p? + {evkf 

< e2-7(^(£)-i I v{vy^\vk\^F{v)dv 

J\v\<M 

< |A:|V-T(^(e)-i [ \v\^iy{v)-^ F{v)dv 

J\v\<M 

< C \k\^ e^-^ ifie)-^ 

(using the assumption (A2) that jvp v{v)^^ F is locally integrable) and 

It is readily seen (using (j46p and the fact that 7 < 2) that 

dl [p, k) — > as e ^ 0. 

Furthermore, d\ is bounded, for e small enough by C[A;p. So we only need 
to evaluate the limit of ^Kp, A;). For that purpose, we first rewrite as 
follows: 

d^i^k) = ^l {\vr^evkf V{v)F,{v)m) 

^2{P,>^) e-Zifie) Ju>M mv) + \v\-^e"'ip{e)py + i\v\-^evk)^ \v\^+'^-^ 
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where 

u{v) = Iv]-'^ iy{v) and Fo{v) = \v\^^" Fo{v). 

Note that Assumptions (B1-B2-B3) (see also ()28p ) imply that ? and Fq are 
uniformly bounded from above and below for |u| > M for a suitable M > 0, 
and that 

lim ^{v) = uq and lim Fq(v) = kq. (44) 

t;|— >oo \v\^oo 

We now do the change of variable 

w = e\k\\v\~^v, 

for which we have: 

\v\ = ^— V = —, and dv = ^ dw. 

{e\k\)^-fi, |tt;|i-'5(e|A:|)i-/3 {e\k\)^-i^ \w\ ^-f^ 

We obtain (with e = k/\k\): 
dl{p,k) = 



J\^\>M^-Pe\k\ {y^{w) + \w\'^P{e\k\)^e^^{e)pY + {wef 



\w\ 1-f V(el^l)'-''/ 



where 

w 



and 



-13 1 

1-/3 (e|A;|) 1-/3 



V|t«| 1-/3 (e|/cj) 1-/3 / 

(we dropped the dependence in k in and -Fq to keep notation simpler). 
The definition of 7 thus gives: 

dl{p,k) = 

[w ■ e)^ 



{l-P)-^\k\^ [ ^ ^ 

i|«)|>Afi-'3e|A:| (i/e(it;) + i-/3 |A:| 1-/3 y?(e)e i-/3p)2 + (tu • e) 



2 



^I7^HF|(^) / Hi-. 1 
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Finally, we note that for any w and k, we have (note that /? < 1) 

lim D^{w) = vq 

£—*0 

lim Fq (w) = kq 

£—>0 

and 

lim ( = hm -^i ( = 1. 

Thus, the integrand in (jl5|) converges pointwise to 

(w ■ e)^ z/pAtQ 
i/2 + (w;.e)2 |u;|^+7- 

Furthermore, it is bounded by 

i/f + (u- . e)2 1^ ^ V |fc| 

thanks to the following lemma and an appropriate choice of M, see the 
following lemma (extracted from [3l Theorem 1.5.6]): 

Lemma 4.4 (Potter's Theorem). For all 6 > there exists some constants 
C, M > such that 

< C(l + lAI"^) for alls>M and A > M/s. 

Moreover for any C > 0; 

/£(s) — > oo and s-^^(s) — > as s — > oo. (46) 

Choosing < (5 < 7/2 (and thus fixing a corresponding constant M), we 
deduce that 

l4(p,A;)| <C7|A:P(1 + |A;|-^) 
and by Lebesgue's dominated convergence theorem we get: 

Jm.n u^ + {w eY \wY^+^ 

which concludes the proof.Q 




l{Xs) 
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4.4 The remaining term 



We will also need the following technical lemma to pass to the limit in the 
first and last term in (|39p . Note that we directly include in this lemma the 
limit case 7 = 2 since we shall need it this same technical lemma in this 
limit case for the proof of Theorem 13.41 

Lemma 4.5. Assume (A1-A2) and (B1-B2-B3) with a > 0, /3 < a and 
P <2-a. 
Then 



V + Qie)]) + eiv ■ k 



1 



iy{v)F{v) dv < C{p, k) e>^W7/2+^ ; 1} ^47^ 



for some rj > 0, all k eM.'^ , p e M+, and C{p, k) G L\ 



00 ( 



)■ 



Remarks 4.6. Note that the assumptions of Lemma \4-5\ imply thatO < 7 < 
2, and therefore j = 2 is included in these assumptions, except when a = 1. 
In other words we exclude the case a = (3 = \ and hence when j = 2 we 
always suppose a > 1. 

Proof of Lemma 14.51 

To prove (j47p . we split the integral into three part parts: 



u + 9{e)p + eiv ■ k 



1 



where 



c±ip,k) :-- 



,|±i< 



u{v)F{v) dv = c\{p, k) + ci(p, k) 
v{v) 



V + Q{e) p + eiv ■ k 



ij(v)F(v) dv 



(note that this splitting is well-defined since /? 7^ 1 in the assumptions). 
In order to get a bound for c^_, we note that (|42p implies: 



< d{e) \p\ + e'^u^'^iv -kf + e \v ■ k\ 



V + Q{e)p + eiv ■ k 
and thus 



c+(j3,/c) < 6'(e)|pl J F{v)dv 
+{^\k\f 



-1 ^ Ivl"^ u~'^{v)F{v)dv 



+e|A:| 



_i . \v\ F{v) dv. 

\v\< 
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We write 

c%{p,k) < e{£)\p\ + {e\k\fc%^2iP,k)+e\k\c'+^3{p,k). 

The computation of c^2(P)^) requires the computation of the fohowing 
integral: 

/r -1 . \v\'^u-'^{v)F{v)dv 

Ji \v\< \ 

which we estimate using the assumption (A2) that \v\'^ ^{v)^^ F is locahy 
integrable, and the fact that F{v) < Cii{v) \v\~'^~°' and v{v) > C2 \v\^ for 
\v\ > M. When /3 / 2 - a, it yields 

for some S as small as we want, using the Potter's lemma to estimate i. 
In the case when P = 2 — a, it yields 

4,2<C(1 + H(e-V(i-/^))) 

with 

/■^ i(r) 

E{z) := -:-Ldr if z> M, E(z) := if z < M. 
Jm r 

In both cases this yields, by fixing 6 = 7/2, using the trivial estimate 

z-^E{z) > 

2— > + 00 

for any 5 > 0, and keeping only higher order terms, we obtain the following 
(non optimal) bound: 

{e\k\fc%^2{P,k) < 

for some ?? > 0. 

Similarly we found 

4,3 <C (1 + e^-i-^) if a / 1, 4 3 <C (l + E (e-^/^^-'^))) if a = 1 

for some 5 as small as we want. 

Using the same arguments as before (and recalling that when a = 1, 
7 = 2 is not allowed since (3 < a), we get 

iE\k\)c%^,ip,k)<Cie\k\r^^^/^+^-''^ 
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for some 77 > 0. 

Finally, in order to get a bound on dL, we do the change of variable 
1 

w = £1-/3 V which leads to 



ct{p,k) < L _i . 2v{v)F{v)dv 

\v\>e^ \ 



a-0 



,\>i \w 



N+a-P 



< Ce^ifie) I \w\ 

'\w\>l 



-N-a+l3 



1 + \wY' dw 



where we used Potter's lemma for some (5 > 0. 

Collecting all the terms, and keeping only higher order terms, we find 
easily (I17|. □ 

Finally we show that the last term in ([39|) goes to zero in P'([0, 00] x M^): 
First, we note that for any g S Lp'{vF~^), we have by Cauchy-Schwarz 
inequality and assumption ([27j) 



K{g)\ < / a{v,v'Mv')\dv' 



< 



F' 



1/2 



^ dv' 
< M u F \\g\\L2^^p-iy 
Next, using Lemma [^3} we deduce that 



V dv 



F{v') 



1/2 



1 



9{e) 



< 



C 



i^{v) + 9{e)p + eiv ■ k 
u{v) 



K{g^){v) dv 



iy{v) + 0{e)p + eiv ■ k 



u{v)F{v)dv\ Wg'^WL^uF-i) 



C{p, k) u-e\\ 



Note that here since 7 < 2, one can find 77 > such that 7/2 + ?] < 1, and 
therefore the estimate ()47p in Lemma 14.51 yields a bound from the below of 
the form Ce"</'^+'^. 
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Last, recalling inequality ([38|) . we have for any test function cj) G I'((0, oo) x 



/o JR^ Jm^ 



v{v) + ^(e)p + eif • k 



K{g^){v) dv dk dp 



< 



c{p,ky 4>{p,ky dk dp 



1/2 



< 



^ ll5'^llL°°({a,oo);L2(i^F-i)) 

C{4>) 



9 IIl2{(0,oo);L2{!/F-1)) 



which goes to zero thanks to (06]) and the fact that r/ > 0. □ 



4.5 Conclusion 

We now have all the tools to pass to the hmit e ^ in (j39p . The first part of 
Lemma 14.51 imphes that the first term converges to po{k) locally uniformly 
with respect to k and p, while the second term can be rewritten as a^(p, k)p^ 
which converges, using Lemma 14.31 in T>'{{0, oo) x R^) to {—p — K\k\'^)p. 

The remaining term goes to thanks the the second part of Lemma 14.51 
and the preceeding subsection. 

Hence we can pass to the limit in (I39p and recover (|2ip . The proof of 
Theorem 13.21 can then be completed as the proof of Theorem 11. 1[ 



5 Proof of Theorem 13.41 

The proof of Theorem 13.61 is very similar to the proof of Theorem 13.21 We 
recall the main steps for the reader's sake. 

We recall that we define the following time scale in the critical case: 

e{e) = £^£{£'^) In(e-i). 
and we shall use the same notation as in the previous section: 
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The starting point is again the equaUty (j39p . As in the previous section, 
we see that if we take that the first term converges to po when e goes to 
zero. 

Proceeding as in Section 14.41 we can also show (using Lemma 14.51 with 
7 = 2) that the last term goes to zero: 



1 



e{e) 



v{v) 



< 



c 



v{v) + 9{e)p + eiv ■ k 



1 



v{v) + 0{e)p + eiv ■ k 



v{v)F(v) dv 



\9 \\l'^(uF-^) 



Note that here since 7 = 2, the estimate (j47p in Lemma 14.51 yields a bound 
from the below of the form C e. 

Last, recalling inequality ([38|) . we have for any test function G I'((0, oo) x 



Jm.^ ^{^) JR^ 



u{v) 



1 



K{g^){v) dv dk dp 



< 



< 



Ce 

W) ^ 

9{e) 



R^ 



u{v) + 9{e)p + eiv ■ k 

C{p,kf (t){p,kf dk \ dp \ ||ff^||Lo°((a,oo);L2(i.F-l)) 



\9 IIl2((o,oo);L2(j.f-1)) 



< 



C{<t^) 



^{e) ln(e~i) 

which goes to zero thanks to assumption (j32p in Theorem 13. 4t l{r) Inr ^ 
+00 as r ^ +00 indeed immediately implies that £(e~^/(-'^~^)) ln(e^-'^) 
+00 as e ^ (recall that here 1 — f3 = a — 1 >0). 

Finally, we are left with the task of determining the limit of the symbol 

in 

It has already been proved that the first term in the right hand side of 
()43p is bounded (uniformly in e) by \p\ and converges to —pJ^N -^(^) dv = 
—p as e goes to zero. 

So it only remains to show that 

{evkf 



d%p,k) :-- 



1 



{u + 9{e)pY + {ev ■ k)' 



v{v) F{v) dv 



converges to k |A;p and is locally bounded when e goes to zero. 
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For some M > 0, we write 

d'{p,k) = dl{p,k) + d'2{p,k), 



where 



d\(jp,k) = —— / — — —-rv[v)F{v)dv 

e{e) J\^\<M{^{v) + G{e)pY + {evkY 



^ / M / -n / v{v)-^\v-k\^F{v)dv 



< 



< 



C\k\ 



2 



(using the assumption (A2) that \v^ v{v)^'^ F is locally integrable) and 

dl{p,k) = ^ f ' l^' t,(^y)F{v)dv. 

^(e) J\v\>M {'^{v) + 0{e)py + {ev- ky 



Concerning (p, k), as above we deduce from the assumption ([32]) on I 
that 

ip{e) ln(e^"^) +oo as e ^ 0, 

and therefore df (p, k) goes to as e goes to 0. Moreover it is clearly bounded, 
for e small enough, by some C 

Now it remains to evaluate the limit of did?, fc). For that purpose, we 
first rewrite d^ as follows: 

dHv k)-^ I {\v\-^evkf 

"^2^^' - e{e) y|„|>M {V{v) + \v\~P 9{e) pf + {\v\~P evkf \v\^+--P '^'''1^ 
where 

V{v) = \v\-l^ v{v) and Fq{v) = Fo{v). 

Note that Assumptions (B1-B2-B3) (see also (|28p ) imply that v and Fq are 
uniformly bounded from above and below for |?;| > M and that 

lim ^{v) = uq and lim Fo{v) = kq. 

\v\^oo \v\^oo 
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We do again the change of variable w = £\k\\v\ (recall again that (3^1 
from the assumptions), and we obtain (with e = k/\k\): 

dl{p,k) = 

(l-/3)-i [ {wef 



^ ^+LV (g|fe|) ^ ' ' ^ (e|^l) H '-^ dw, 
\w\^^ \{e\k\)^-fi J 



where 

w 



and 



-13 1 

|'w|i-'3(e|A:|)i-'3 



^0 H = ^0 — I 

V|w|i-'3(e|A:|)i-/3 / 

(we dropped the dependence in k in and Fq to keep notation simpler) 
The fact that 7 = 2 now yields 

dl{p,k) = 

(1-/3)"' ,,,2 f {w-e? 



^i7^(.)Fo^(.),^ (43) 



\w\N+2 



\{e\k\) 
Lemma |5. II below implies 



d\{p,k) K\k{' 



which concludes the proof of Theorem [37 

Lemma 5.1. We have: 

(l-/3)-i /• {wef 



^{e) In(e-i) y|^|>Mi-.,|,| j^-.^^^ ^ |^|_i||^|_^^(^)^_^^^j2 ^ _ 



e)2 



1 



?nz.)Fo^(.) J H-^ (49) 



1-13 



as e ^ 0. 
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Proof of Lemma I5.lt First, we show that 



is well defined: We denote 



and then we have 



|^„|=A '"O+Wf \W\ 

' aa[w) 



XN+2 



|«,|=A 1^0 + ^1 



da{v) 



We deduce from L'Hopital's rule that 



' Hm^ = ^lim^= / 



(1 - /?) A^o In(A-i) (1 - (3) A^o In'(A-i) (1 - (3) i^o J\v\ 



=1 



Then, we note that from the assumptions (B1-B2) and the definition of 
a slowly varying function, for any w and k, we have 

u^{w) = f + o(l) as e ^ 0, 

Fq{w) = kq + o{l) as e ^ 0, 



and 



e ' ^ , = ' ' 1 = 1 + o(l) as e ^ 0. 



Thus, the integrand in ([^8]) divided by (p{e) converges pointwise to 
{w ■ ef uqKo 



u^ + {w- e)2 |u;|^+2 
and it is bounded by 

^ {^-e? 1 



as e — > 0, 



N+2 • 
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Then, easy but tedious computations yield, for bounded p and k and for 
e going to 0: 

{w ■ e)^ 



\w\>M^-0e\k\ ^ \w\^\k\~^{e)e^pj + {w ■ e)2 



- / ~^ 2l — uvT^dw 

J\w\>iin-Pe\k\ H^w{ |u;|^+2 

7|«,|>Mi-/3£|fe| 1^0+^1 |t(;r+2+w(i W 

< o(l) In(e-i) + C(/7(e)e2 

< o(l) ln(e"^). 

Thus we have 

1 /■ («; • e)2 



hm 



|y;|7V+2 



dw 

.{e\k\)^J 

,,2 



hm — / — H K — m-nr dw. 



Finally, using the fact that for all k ^ 0, 

, ln((M^-f^e\k\)-^) 

hm — — = 1, 

we deduce ()49p , which completes the proof of Lemma 15.11 □ 

6 Appendix 

Lemma 6.1. Let us assume that b = b{v, v') = \v — v' f with P G 
(i) For (3 G (— A^, a) there exists C such that 

{v)^ < u{v) := K{F){v) < C {vf, ^vGR^. 
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(ii) For (3 € (— min{a; A^}, min{Q; A^}) there exists M such that 

F'^dv'<M, VueM^. 
5 

(in) For (3 G (— min{a; A^/2}, a) there exists M such that 

pi ip. 

— ^dv'<M VvGM^. 



Proof of Lemma 16.11 Point (i). The case (3 > being very simple, we 
only deal with the case /3 < and large values of |f |. We split v = vi + 1^2 
with 

vi{v)= f F{v')\v-v'f dv' , i^2iv)= [ F{v')\v-v'f dv'. 

J\v'-v\<\v\/2 J\v'-v\>\v\/2 

On the one hand 

ui{v) < F(\v\/2) [ \v-v'fdv'<Fo\v\-''-^\v\^+'^<C\vf. 

J\v'-v\<\v\/2 

On the other hand 

'^2{v) < {\v\/2f [ F'dv' <C\vf. 

J\v'-v\>\v\/2 

Finally, 

J^2{v) > {S\v\/2f f F'dv' > {3\v\/2y'^ [ F'dv' > C'^ \vf 

for \v\ > 2. 

Points (ii) and (Hi). From Lemma l6. II we have on the one hand 

^dv'=[ F'\v-v'\-^ <C{v)-l^ <Mv-'^ V^;eM^, 



and on the other hand 



^b^dv'< I {v')-''-"-('\v-v'f'' <C{vff' <Mu\ Vt;€™^ 



□ 
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